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ABSTRACT 

Using an XMM-Newton observation of the Perseus cluster we show that the 
excess in the flux of the 7-8 keV line complex previously detected by ASCA and 
BeppoSAX is due to an overabundance of Nickel rather than to an anomalously 
high Fe He/3/Fe Hea ratio. This observational fact leads to the main result 
that resonant scattering, which was assumed to be responsible for the supposed 
anomalous Fe He/?/Fe Hea ratio, is no longer required. The absence of resonant 
scattering points towards the presence of significant gas motions (either turbulent 
or laminar) in the core of the Perseus cluster. 

Subject headings: X-rays: galaxies — Galaxies: clusters: individual: Perseus 
Abell 426 — Galaxies: abundances — intergalactic medium — scattering 
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1. Introduction 

The X-ray emission from clusters is due to a diffuse, tenuous (with typical densities 
of 10~^ — f0~^cm~3), and hot (with typical temperatures of 10'' — 10^ K) thermal plasma. 
Although for these ranges of density and temperature the gas is optically thin to Thomson 
scattering for the continuum, it can be optically thick in the resonance X-ray lines of highly 
ionized atoms of heavy elements (Gilfanov et al. 1987). Apart from other interesting observ- 
able effects (Sazonov et al. 2002), the major effect of resonance scattering (the absorption of 
a line photon followed by immediate re-emission) is to distort the surface-brightness profile 
of the cluster in the resonance line due to diffusion of photons from the dense core into the 
outer regions of the cluster. This must be taken into account when attempting to determine 
element abundances from X-ray spectroscopic observations of galaxy clusters. In fact only 
with the key assumption that the plasma is optically thin lines equivalent widths can un- 
ambiguously convert to element abundances, when fitting CCD spectra with the available 
plasma codes. In presence of resonant scattering the true abundances in the core of clusters 
are significantly underestimated because the line emission is attenuated due to photons scat- 
tered out of the line of sight. To make things worse, the most promising hne for resonant 
scattering is the Hea Fe emission line at 6.7 keV (Gilfanov ct al. 1987; Sazonov et al. 2002) 
which is also one of the most prominent emission line in cluster spectra and in general drives 
the global abundance determination. 

High sensitivity and high resolution spectrometers are needed to directly measure the 
spectral features of resonant scattering (as modification of the line profile or resolution of 
the He-like line into its constituents in order to determine directly the effects of scattering) 
and polarimeters, to detect the polarized scattered radiation (Costa et al. 2001; Sazonov et 
al. 2002). Currently the simplest method to reveal and estimate the presence of resonance 
scattering is to compare the fluxes of an expected optically thick line and of an optically 
thin one and to check if it is correctly modeled by a plasma code assuming optically thin 
emission. This was done in the past with ASCA and BeppoSAX for the ratio between Hea 
Fe line at 6.7 keV and the He/3 Fe hne at 7.90 keV (which is expected to have an optical 
depth typically smaller than one for resonant scattering) and in particular the best data were 
the ones for the Perseus cluster. 

Molendi et al. (1998) analyzed data collected with the MEGS on board BeppoSAX 
and found that the ratio of the flux of the 7-8 keV line complex to the 6.7 keV hne was 
signiflcantly larger than predicted by optically thin plasma code and that the ratio decreases 
with increasing cluster radius. They noted that this effect could be explained either by 
resonant scattering or by a Ni overabundance, eventually favoring the former explanation. 
On the contrary Dupke & Arnaud (2001), according to the experimental evidence of a central 
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enhancement of SNIa ejecta in cD clusters, favored the over- abundant Ni explanation. 

These were the two hypothesis that the resolution and sensitivity of past instruments 
could not resolve. XMM-Newton has now for the first time the combination of resolution 
and effective area at high energies to give an unambiguous answer to the question. In this 
paper our aim is to try to solve this controversy. A complete spectral and spatial analysis, 
in particular of the temperature structure of the Perseus cluster which requires a detailed 
temperature map, is beyond the scope of this paper. 

The outline of the paper is as follows. In section 2 we give information about the 
XMM-Newton observation and data preparation. In section 3 we present spatially resolved 
measurements of temperature and Ni and Fe abundances. In section 4 we discuss our results 
and draw our conclusions. 

At the nominal redshift of Perseus (z=0.0183), 1 arcmin corresponds to 22.2 kpc {Hq — 
70 km Mpc~^, fl^a = 1 — = 0.3). In the following analysis, all the quoted errors are 
at la (68.3 per cent level of confidence) unless stated otherwise. 

2. Observation and Data Preparation 

The Perseus cluster was observed with XMM-Newton (Jansen et al. 2001) during Revo- 
lution 210, with the THINl filter and in Full Frame Mode, for 53.6 ks for MOS and 51.2 ks 

for PN, but resulting in an effective exposure time (as written in the keyword LIVETIME of 
the fits event file) of 53.1 ks for the MOS and 24.7 ks for the PN. We generated calibrated 
event files using the publicly available SASv5.3.3. 

To fully exploit the excellent EPIC data from extended and low surface brightness 
objects and from this observation in particular, the EPIC background needs to be correctly 
taken into account. The EPIC background can be divided into a cosmic X-ray background, 
dominant below 2-3 keV, and an instrumental background, dominant for energies higher than 
2-3 kcV (for what concern the continuum emission, apart from fluorescence lines). This latter 
component can be further divided into a detector noise component, present in the low energy 
range (below 300 eV) and in a particle induced background, which is the major concern 
for our scientific case. The particle induced background consists of a flaring component, 
characterized by strong and rapid variability, produced by soft protons (with energies less 
than few hundreds of keV) which are funneled towards the detectors by the mirrors, and 
a second more stable component associated with high energy particles interacting with the 
structure surrounding the detectors and the detectors themselves. The latter component has 
been studied using CLOSED filter observation, shows only small intensity variations and is 
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characterized by a flat spectrum and a number of fluorescent lines. Apart from a rather 
strong variability of the fluorescent lines this component can be properly subtracted using 
a large collection of background data. The common way to face the flaring component is 
to remove periods of high background, because the S/N is highly degraded, especially at 
high energy (where the data are crucial to measure the exponential cut-off and thus the 
temperature of the emitting plasma) and because the shape of the spectrum is varing with 
time (Arnaud et al. 2001). The strategies to reject these flaring periods are mainly two: 
selection of time intervals where the count rate in a given high energy band is lower than 
a given threshold (which has been our approach in previous analysis where we flxed the 
thresholds at 0.35 cts/s for PN in the 10-13 keV band and 0.15 cts/s in the 10-12 keV band 
for MOS, based on Lockman Hole data) or flnding a mean count c and then choosing as a 
threshold value c + Scr, by means of a Gaussian or Poissonian fitting or o"-clipping methods 
(see Appendix A of Pratt & Arnaud 2002 for the second approach and Marty et al. 2002 for 
a general discussion on soft proton cleaning criteria). 

The light curve in the 10-13 kcV band for the PN observation of the Perseus cluster 
is shown in Fig.l together with our standard threshold of 0.35 cts/s. It is evident that the 
observation is badly affected by soft proton and if we adopt our threshold all the observation 
would be rejected. The light curve is also structured in such a way that a 3(j-clipping method 
rejects only 529 s of observation, flnding a mean rate of 0.53 cts/s with a standard deviation 
of 0.15 cts/s, while fltting with a Gaussian and rejecting all the intervals above 3a from 
the mean rejects only 700 s of observation. Our approach was therefore to consider all the 
observation for two reasons: we can exploit the fact that Perseus is the brightest X-ray 
cluster and it is so bright in its central zone that the background, also in presence of a high 
level of soft protons as we have in our observation, is not important; moreover we can try to 
model the soft proton which contaminate the spectra using in first approximation a power 
law as a background model (which means that the model is not convolved via the effective 
area of the instrument). The self-consistency and viability of our approach will be shown in 
the results. 

We have accumulated spectra in 9 concentric annular regions centered on the emission 
peak with bounding radii 0.5' - 1', 1' - 2', 2' - 3', 3' - 4', 4' - 5', 5' - 6', 6' - 8', 8' - 10', 
10' — 14'. We did not consider the inner bin inside 0.5' in order to avoid contamination by 
the power law spectrum of the Seyfert cD galaxy NGG 1275. 

Spectra have been accumulated for the three cameras independently and the blank fields pro- 
vided by the calibration teams were used as background (Lumb 2002). Background spectra 
have been accumulated from the same sky regions as the source spectra, after reprojection 
onto the sky attitude of the source (this ensures the proper subtraction in the same way as 
it was performed in detector co-ordinates, see Lumb 2002). 
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The vignetting correction has been apphed to the effective area generating effective area files 
for the different annular regions using the SAS task arfgen. We generate flux weighted arf 
using exposure corrected images of the source as detector maps and the parameter extended 
source switched to true, following the prescription of Saxton & Siddiqui (2002). Spectral 
results for the cluster A3528 obtained in this way and with the vignetting correction ap- 
plied directly to the spectra (Arnaud et al. 2001) are practically the same (Gastaldello 
et al. 2003). We also correct the PN spectra for out of time events following the pre- 
scriptions of Grupe (2001). The redistribution matrices used are ml_r6_all_15.rmf (MOSl), 
m2_r6_alL15.rmf (M0S2) and, depending on the mean "RAWY" of the region, the set of ten 
single-pixel matrices, from epn_ff20_sY0.rmf to epn_ff20_sY9.rmf, and double-pixel matrices, 
from epn jf20_dY0.rmf to epn_ff20_dY9.rmf, for PN. 

Due to its higher effective area (further increased by the use of doubles data) and similar 
spectral resolution at high energies, the PN camera will be the leading instrument in our 
analysis and the one for which the results are most compelling, in particular for what con- 
cerns the Ni abundance. There are still some problems for what concern the three EPIC 
cameras cross-calibration and in particular at high energies the study of power-law sources 
returns harder spectra for MOSl, intermediate for M0S2 and then the softest for PN (Kirsch 
et al. 2002). Also our analysis of the galaxy cluster A3528 gives systematically higher tem- 
peratures and abundances for MOSl respect to M0S2 and PN. The conclusions of a recent 
work aimed at assessing the EPIC spectral calibration using a simultaneous XMM-Newton 
and BeppoSAX observation of 3C273 strengthen this fact: the MOS-PN cross calibration 
has been achieved to the available statistical level except for the MOSl in the 3-10 keV band 
which returns flatter spectral slope (Molendi & Sembay 2003). 



3. Spectral modeling and energy ranges used 

All spectral fitting has been performed using version 11.2.0 of the XSPEC package 
(Arnaud 1996). 

As a first step we concentrate on the hard band which is the one of interest to de- 
termine the abundances of iron and nickel and also to make a direct comparison with the 
MECS results. We use three different energy bands: 3-10 keV, 3-7 keV in order to have a 
band less contaminated by the hard tail of soft protons, and the 3-13.5 keV and 3-12 keV for 
PN and MOS respectively in order to have more data to acceptably model the soft protons 
background. 

When fitting the first two bands we analyze the spectra with a single temperature VMEKAL 
model (Mewe et al. 1985; Kaastra 1992; Liedahl et al. 1995) with the multiphcative compo- 
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nent WABS to account for the Galactic absorption fixed at the value of 0.143 x lO^^cm"^ 
(according to Schmidt et al. 2002). We leave the abundances of Ar, Ca, Fe and Ni (the only 
elements which have emission lines in the range 3-10 keV) free and keep all the other abun- 
dances fixed to half the solar values (Fukazawa et al. 2000) (this corresponds for example to 
the IT (3-10 keV) model in Tab. 1 and Tab. 2). 

When fitting the wider high energy band, more contaminated by soft protons, we add a 
power law background model (VMEKAL-^POW/B in XSPEC) in order to model the soft 
proton background component (this corresponds to the lT-|-pow/b (3-13.5 keV) model in 
Tab. 1 and IT+pow/b (3-12 keV) Tab. 2). 

As a second step we fit the entire energy band 0.5-10 keV with two models: 
a single temperature model leaving Nh to vary freely (the fit is substantially improved respect 
to the one with Nh fixed to the galactic value) and the abundance of O, Ne, Mg, Si, S, Ar, 
Ca, Fe and Ni. For the outer annuli, when required from the previous analysis in the hard 
band, we add the pow/b component with normalization and slope fixed at the best fit values 
found (these models corresponds to IT (0.5-10 keV) or lT-|-pow/b (0.5-10 keV) in Tab. 1 
and Tab. 2); 

a two temperature model ( WABS*(VMEKAL+VMEKAL) in XSPEC) where the metal 
abundance of each element of the second thermal component is bound to be equal to the 
same parameter of the first thermal component. As for the single temperature model we 
add the pow/b component when required (these models corresponds to 2T (0.5-10 keV) or 
2T+pow/b (0.5-10 keV) in Tab. 1 and Tab. 2). 

The two temperature model is a rough attempt to reproduce the complex spectrum resulting 
from projection effects, azimuthal mean of very different emission regions (like holes and 
luminous regions in the Perseus cluster, see Schmidt et al. 2002; Fabian et al. 2002) and an 
atmosphere probably containing components at different temperatures, as in M87 (Kaiser 
2003; Molendi 2002). 

We also allow the rcdshift to be a free parameter in order to account for any residual 
gain calibration problem. We adopt for the solar abundances the values of Crevesse & Sauval 
(1998), where Fe/H is 3.16 x 10"^ To make comparison with previous measurements, a 
simple rescaling can be made to obtain the values with the set of abundances of Anders & 
Crevesse (1989), where the solar Fe abundance relative to H is 4.68 x 10~^ by number. 
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4. Results 

4.1. IT results in the high energy band 

In Fig. 2 we show the temperature profile obtained analyzing the single events spectrum 
for the PN camera. This is also an example of our working procedure. The full circles 
refer to the results obtained using the 3-10 keV band, while the open triangles indicates 
the results obtained using the 3-7 keV band with the Ni abundance frozen to the best fit 
value obtained in the 3-10 keV band. It is clear that where the source is overwhelmingly 
bright the hard component of the soft proton does not affect the spectrum and there are 
no differences between the temperatures obtained in different energy bands, while in the 
outskirts of the cluster, where the source brightness is lower and the soft protons become 
important, the fitted plasma temperature reaches uncorrect and unphysically high values and 
large residuals at high energy are present. With the open squares we show the temperature 
obtained by fitting not only the source but also the soft protons with a power law background 
model, in the energy band 3-13.5 keV: as we expect in the inner region adding the background 
component does not affect the temperature determination nor statistically improve the fit, 
on the contrary in the outer regions the temperature are significantly reduced and the fit 
is improved, eliminating the residuals at high energies. For example in the 6' — 8' ring the 
simple single temperature fit gives a of 1502 for 1137 degrees of freedom, while the fit 
with the power law background model in addition gives a of 1258 for 1288 d.o.f. 

To confirm our results we compare the temperatures obtained in this way with those 
obtained with the MEGS instrument on board BeppoSAX (De Grandi & Molendi 2002a). 
The temperature profiles, apart from the differences in the three camera due to the cross- 
cahbration problems we discuss before (confirmed also with the superb statistics of Perseus) , 
are in good agreement at least up to 8 arcmin. In the outer rings between 8 and 14 arcmin 
the increasing importance of background relative to source counts prevent us from recovering 
a correct temperature with our method (see De Grandi & Molendi 2002a for a more general 
discussion about XMM-Newton and BeppoSAX temperature determinations and the greater 
sensitivity of the latter over the former to low surface brightness regions due to much lower 
background) . 

With a determination of the temperature structure we can address the issue of metal 
abundances measure and attempt to discriminate between the presence of resonant scattering 
or the supersolar abundance of Nickel. Resonant scattering is increasingly important towards 
the center of the cluster so we choose our two inner bins to test its presence. Fitting the 
spectra with a MEKAL model, assuming solar ratios, actually does not reproduce the 8 
keV line complex. As shown in the left panel of Fig. 4 for the 1' — 2' bin, the emission is 
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underestimated as for previous missions (see Fig.l of Molendi et al. 1998 for example). 
However the data show for the first time that the excess is due to an uncorrect modehng of 
the Ni Hea hne complex at 7.75-7.80 keV (in the rest frame of the source) and not to an 
underestimation of the Fe B.e(3 line which is correctly modeled. Infact if we fit the data with 
a VMEKAL model, we eliminate almost completely the residuals and give a better fit with 
a Ni abundance of 1.23 in solar units, as shown in the right panel of Fig.4. The fit with a 
MEKAL model gives a of 855 for 802 d.o.f for the 0.5' - 1' bin and 1116 for 1023 d.o.f. 
for the 1' - 2' bin, while a fit with a VMEKAL model (with Ar and Ca fixed to 0.5 Z/Zq, 
because they are not important in driving the fit, in order to have only the Ni abundance 
as additional free parameter) gives a of 835 for 801 d.o.f. for the first bin and 1092 for 
1022 d.o.f. for the second bin, with Ax^ which arc statistically significant at more that the 
99.9% according to the F-test (the value of the F statistics is F=19.2 with a probability of 
exceeding F of 1.4 x 10^^ for the first bin and F=22.5 with a probability of exceeding F of 
2 X 10~^ for the second bin). 

We can conclude that the ratio of He/S/Heo; Fe lines is not anomalously high respect to 
the optically thin model and that it is not necessary to invoke resonant scattering in the core 
of the Perseus cluster. The excess in the flux in the 8 keV line complex respect to a MEKAL 
model is entirely due to Ni overabundance with respect to solar values, as was previously 
suggested (Dupke & Arnaud 2001). 

The reader will notice some residuals in the Hea Fe line complex at 6.7 kcV. This is an 
instrumental artifact present only in the inner bins out to 2' of the FN camera wc suspect 
connected to some residual CTI problems due to the high flux of the Perseus cluster. The 
net effect is to lower the energy resolution broadening the fine profile. We test that this does 
not affect our results fitting spectra for our two inner bins with a bremsstrahlung model 
plus two Gaussians fixed at the energies of the Fe Hect at 6.67 keV and Fe He/3 at 7.90 keV 
leaving the redshift, width and normalizations of the two lines as free parameters. We find 
that the Gaussian width of the Hea line in the two bins are 5.3 x 10~^ keV in the 0.5' — 1' bin 
and 4.5 x 10^^ keV in the 1' — 2' bin. If we force the He/5 to have up to a width of 8 x 10^^ 
the excess due to the Ni He« line blend at 7.75-7.80 keV is still significantly present. This 
instrumental effect is evident because of the large equivalent width of the Fe line at 6.67 keV 
and does not alter significantly the measure of metal abundances as we show further on. 

We can make some other important considerations investigating another line ratio, 
namely the Hca Fe line complex at 6.7 kcV over the Ha Fe line at 6.97 kcV. This ratio 
allows a robust and independent determination of the temperature, because as the temper- 
ature increases the contribution from the He Fe fine decreases while the contribution from 
the H Fe line increases. Thus the intensity ratio of the two fines can be used to estimate 
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the temperature. This was done in the past determining the variation with the temper- 
ature of the centroid of the blend of the two hnes, because gas proportional counters did 
not have sufficient spectral resolution to resolve the two lines (Molendi et al. 1999). Now 
with XMM-Newton we can resolve the lines, measure separately their intensity and use their 
ratio as a thermometer. To do that we obtain a calibration curve of the line flux ratio as a 
function of temperature simulating spectra with MEKAL model and the PN singles response 
matrix with a step size of 0.1 keV, fixing the metal abundance of 0.3 solar units and the 
normalization to unity in XSPEC units (however the flux ratio is independent from these 
two quantities), with an exposure time of 100 ks to ensure neghgible statistical errors. We 
then model the spectra with a bremsstrahlung model plus two Gaussians for the two iron 
lines, in the energy range 3-10 keV and obtaining the fluxes of the two lines from the best 
fit models. We obtain a calibration curve identical to that of Nevalainen et al. (2003). We 
then measure the line flux ratio from the cluster PN singles data using the energy range 
5.0-7.2 keV to minimize the dependence from the continuum and calibration accuracy and 
to better describe the lines. We fltted each spectrum with a bremsstrahlung model plus two 
Gaussians (using ZBREMSS plus two ZGAUSS models in XSPEC) leaving all the parame- 
ters free, included the redshift (to take into account any possible gain calibration problem), 
except the line energies. The fits for all the annular bins were good with a reduced never 
worse than 1.1 and the results for the temperature derived from line flux ratio are plotted 
as diamonds in Fig. 2. 

As we can see also this independent temperature determination is in good agreement with all 
the others at least out to 3 arcmin where the cluster is very bright and in good agreement out 
to 8 arcmin with the measurement obtained from the model with the power law background 
component, confirming the validity of our modeling. In the last two bins the temperature 
derived from the lines ratio agrees well with the MEGS measurement and starts to differ 
from the determination with power law background component, pointing to the fact that 
our modeling is not sufficient to fully take into account the background in these bins where 
the source is too dim compared to the soft proton background. We can conclude that our 
temperature determination is reliable out to 8 arcmin. 

The concordance between lines ratio and continuum temperature determination adds another 
piece of evidence against resonant scattering. In fact since the Fe Ha line optical depth is 
1.8 times smaller than the Fe Hea one (this is the difference in their oscillator strength), if 
resonant scattering is present, we would expect the ratio of Hea/Ha lines to be lower than 
in the optically thin case. In turn this would lead to an overestimate of the temperature. 
Since this is not the case we can conclude that resonant scattering is not present. 

In Fig. 5 we plot the abundance profiles of Fe and Ni determined by our best fit model 
(thermal model plus power law component for the soft proton background, in the 3-13.5 keV 
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band for PN and 3-12 keV for MOS). We find an evident gradient in both elements: for Fe 
it agrees well with previous determinations, as the BeppoSAX one (without considering the 
corrections for resonant scattering, as done in Molendi et al. 1998), while we have for the 
first time a detailed abundance gradient for Ni, with measurements reliable out to 8 arcmin 
(wc show only the PN data as we discussed before). We stop at this radius because at larger 
radii the temperature determination is no longer reliable and strong emission lines of Ni, 
Cu and Zn induced by particle events affect the spectrum in the crucial range 7.5-8.5 keV 
(Preyberg et al. 2002). It is evident that there are some problems with the iron determination 
by MOSl, as we also found in A3528. 

Knowing that the excess in the 8 keV line complex is due to the Ni line, we can go back 
to BeppoSAX-MECS data and fit them with a VMEKAL model allowing Ni abundance to 
be free. We find an abundance profile in agreement with the more detailed XMM-Newton 
one, as shown in Pig.5. 

4.1.1. Results with the APEC code 

We try to cross-check the results obtained with the MEKAL code with the ones obtained 
with the APEC code (Smith et al. 2001). Churazov (2003) pointed out that the APEC code 
has different energies for the Ni Hea line complex, resulting in a different fitting for the line. 
For comparison with Fig.4 we show the fit of the V — 2' bin with the recently released APEC 
version 1.3.0 in Fig. 6 (but we have to notice that differences in the redshift determination 
could play a role too: the fit with MEKAL found a redshift of 1.49 x 10^^, while the fit with 
APEC found a redshift of 1.54 x 10~^). The excess is still present, although the statistical 
improvement of a variable Ni abundance respect to ratios fixed at solar values is not as 
evident as in the MEKAL case: the fit with an APEC model gives a of 844 for 802 
d.o.f for the 0.5' - V bin and 1138 for 1023 d.o.f. for the 1' - 2' bin, while a fit with a 
VAPEC model gives a of 843 for 801 d.o.f. for the first bin and 1124 for 1022 d.o.f. for 
the second bin (with a probabihty of exceeding the F statistic of 3.8 x 10~^ for the second 
bin). These results are not conclusive, because there are additional issues with APEC, even 
in this latest release: forbidden and inter-combination lines of the He-likc Ni triplet are 
missing and, even after adding these lines, the total Ni emission is still underestimated, thus 
worsening the differences between APEC and MEKAL and making the excess less evident 
(see http : / /cxc . harvard . edu/atomdb/ issues_caveats . html) . 
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4.2. IT and 2T results in the 0.5-10 keV band 

We fit one and two temperature models to M0S2 and PN single data (we avoid MOSl 
data for the calibration problems explained in the previous section) in the full energy band 
0.5-10 keV. Single temperature models cannot adeguatly fit the entire band spectra, giving 
temperatures systematically lower than the ones obtained in the hard band and leaving 
large residuals at high energies. These facts hint towards the presence of more than one 
temperature component, in-fact a two temperature model yields a substantially better fit 
than the one temperature model, although it is still not statistically acceptable, as the 
reduced shows in Tab. 1 and Tab. 2. 

The temperature profiles for PN singles data and M0S2 are shown in Fig. 7: the two 
temperature fit shows the presence of a hot and of a cold component. The temperature of 
the hot component especially in the outer bins matches the temperature determined with the 
fit in the hard band, while the temperature of the cold component is less constrained, being 
about 2 keV in the PN fit and oscillating between 2 and 3 keV in the M0S2 fit. The relative 
normalization of the two components, shown in Fig. 8, shows that the cool component is 
stronger in the center of the cluster, as we expect for cool core clusters. Although there are 
some puzzling results, as in the inner two bins of the M0S2 data where the fitting procedure 
prefers to give more importance to the cool component, and the presence of cool emission 
also in the outer bins where the emission should be negligible (although the cooling radius 
for Perseus is ~ 6 arcmin, Peres et al. 1998). 

The Fe abundance profile, shown in Fig 9 is not substantially changed and in particular 
the abundance gradient is even more evident. Instead for the Ni abundance profile, shown 
in Fig 10, the evidence for a gradient is not present. In fact adding a linear component 
improves the fit (x^ = 1.9 for 5 d.o.f) respect to a constant (x^ = 1.9 for 6 d.o.f) for the PN 
Ni abundances derived by the IT model in the hard band, while the Ni profile derived by 
the 2T model in the entire band is essentially flat (fltting a constant returns a x^ = 3 for 6 
d.o.f and a linear component does not improve the flt, = 2.9 for 6 d.o.f). We caution the 
reader that the 2T modelization is rather complex, because some not completely justified 
assumptions are made, as for example that the abundances of the two components are equal, 
and there is some degeneracy in the contribution of the two X-ray emission components and 
the soft proton power-law background (see for M0S2 and PN the substantial difference in 
the temperature of the cool component ^). Therefore the derived Ni abundance should be 
taken with some caution. Moreover it is very difficult to explain, in presence of a confirmed 



^With the latest release of SAS, version 5.4.1, which revise the quantum efficiency for the MOS and the 
PN, the agreement between the two detectors should be better, in particular in the low energy band. 
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Fe abundance gradient, a flat Ni abundance profile and a Ni/Fe ratio which increases going 
outward. 



5. Discussion 

Our main result can be summarized as follows: there is no need to invoke resonant 
scattering in the Fe Hea line in the Perseus cluster core and the Fe abundance determination 
with optically thin emission models is reliable. 

Resonant scattering should be important in the core of galaxy clusters, this is partic- 
ularly true for the Fe Hea hne in the core of the Perseus cluster (it has an optical depth 
of 3.3 according to Sazonov et al. 2002). The optical depth of a resonance hne depends on 
the characteristic velocities of small scale internal motions, which could seriously diminish 
the depth r (Gilfanov et al. 1987. sec Mathews et al. 2001 for an example of a detailed 
calculation in presence of turbulent motions). The absence of a clear evidence of resonant 
scattering strongly points towards the presence of significant gas motions. In fact, following 
Gilfanov et al. (1987), the optical depth is 

r-r\l + vlJvlX"' (1) 

where is the optical depth at the line center in the absence of turbulence (in spec- 
troscopy the word turbulence is used for all hydrodynamic motions of unknown pattern which 
cause a broadening of the spectral lines. In hydrodynamics turbulence has a much more re- 
stricted meaning), Vturb is the turbulent velocity of the gas and vpe — {'^kT^,/ MpeY^'^ ~ Cg/S.S 
is the thermal speed of the iron ions and Cg denotes the adiabatic sound speed in the ICM. 
Thus the absence of resonant scattering, r < 1, and assuming = 3.3 (Sazonov et al. 
2002), implies gas motions with characteristic velocities greater than 0.36 c^, i.e a Mach 
number M > 0.36. 

Studies of optical line emission in the central regions of ellipticals reveal chaotic gas 
kinematics typically about 0.2-0.4 of the sound speed in the hot gas (Caon et al. 2000) 
and since small, optically visible line-emitting regions at T ~ 10^ K are likely to be strongly 
coupled to the ambient gas, as some models predict (Sparks et al. 1989) and clear correlation 
between Hq;-|-[N II] and X-ray luminosities suggests, the hot gas should share the same 
turbulent velocities. The first clear example of resonant scattering, acting on the 2p — 3d line 
of Fe XVII at 15.0 A (0.83 keV), has been recently found in the giant elliptical galaxy NGC 
4636 (Xu et al. 2002) using the refiection grating on board XMM-Newton and measuring the 
cross dispersion profile of the ratio between an optically thin emission blend, the two 2p — 3s 
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lines of Fe XVII at 17.0-17.1 A (0.73 keV) and the optically thick hne at 15.0 A. Xu et al. 
(2002) found that if an average turbulent velocity dispersion more than 1/10 of the sound 
speed is added to the assumed Maxwellian the model becomes incompatible with the ratio 
of the 17.1 A/15.0 A lines. We remind that the detection of the resonant scattering is only 
in the inner 1', in-fact the phenomenon does not affect spectra extracted within a full-width 
of 2' (Xu et al. 2002). Another eUiptical, NGC 5044, was observed with the RGS (Tamura et 
al. 2003) and no evidence of resonant scattering was found in spectra extracted in the full 2'. 
If also in this case a cross dispersion analysis where to show resonant scattering, these would 
rule out possible associations at least at these inner scales with optically line-emitting gas, 
because NGC 4636 and NGC 5044 are the most striking examples of chaotic gas kinematics 
in the sample of Caon et al. (2000). 

Another possible source of gas motions is the activity of an AGN, which is now thought 
to be widespread in the core of galaxy clusters and strongly related to hot bubbles, for 
which one of the best cases is indeed the Perseus cluster (Fabian et al. 2002). The induced 
motions could be either turbulent or laminar, as suggested by the recent Chandra and optical 
results in Fabian et al. (2003a,b) (see the discussion about the flow causing the horseshoe Ha 
filament and the derived velocity of 700 km s^^ which for a sound speed of about 1 170 km s^^, 
for a temperature of 5 keV, implies M ~ 0.6 or about the sound waves generated by the 
continuous blowing of bubbles) . AGN activity could explain the lack of resonant scattering 
also in the other best candidate M87 (but see also the discussion suggesting caution for these 
interpretation in the analysis of RGS data for M87 of Sakelliou et al. 2002). 
What is becoming progressively clearer is that resonant scattering effects must be small and 
confined on small inner scales. 

The Fe abundance gradient confirms the general picture of an increase of SNIa ejecta 
in the center of relaxed cD clusters. The Ni abundance, because Ni is almost exclusively 
produced by SNIa, and the presence of a gradient also in this element could be a crucial 
confirmation of this general picture (sec Dc Grandi & Molendi 2002b which report measures 
of Fe and Ni for a sample of 22 clusters observed with BeppoSAX and in particular their 
Fig. 6 showing a segregation between relaxed cD clusters and not relaxed clusters, with the 
formers with greater Fe and Ni abundances respect to latters). The Ni abundance gradient 
is evident in the fit in the high energy band and, looking back at the BeppoSAX data, we 
can attribute the excess in the 8 keV line complex to an increased Ni abundance. However 
the complex thermal structure of the gas prevents us from reaching a robust determination 
of the Ni abundance profile. Detailed temperature and abundances maps are required to 
address this issue. 
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Fig. 1. — PN light curve in the energy band 10-13 keV together with our threshold of 0.35 
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Fig. 2. — PN singles temperature profile. Uncertainties are at the 68% level for one interest- 
ing parameter (Ax^ = 1). Pull circles represent the temperature obtained using the range 
3-10 keV, open triangles the temperature obtained using the range 3-7 keV and open squares 
the temperature obtained by using the range 3-13.5 keV and adding to the source model a 
power law background component. In the last two bins we do not show the temperatures 
obtained in the 3-10 keV and 3-7 keV bands because they are larger than 10 keV. Diamonds 
represent the temperature obtained by the ratio of the fluxes of Hea to Ha Pe lines. 
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Fig. 3. — Temperature profiles obtained with the various EPIC cameras compared with the 
temperature profile obtained with the MECS instrument on board BeppoSAX . Uncertainties 
are at the 68% level for one interesting parameter (A^;^ ~ 1). Crosses represent tempera- 
tures obtained with MOSl, open triangles with M0S2, open squares with PN singles, open 
diamonds with PN doubles and full stars with the MECS. 
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Fig. 4. — PN singles data for the 1' — 2' bin in the 6-8 keV band and the corresponding fit 
with a MEKAL model, on the left, and with a VMEKAL model, on the right, together with 
the corresponding ratios of data respect to models. 



-22- 



EPIC 
0.8 



lECS Fe Abundance profile 



c 



O 



O 
C 
D 

C 
< 



C 

_o 

o 
m 

CD 

u 
c 
o 

T3 

c 

Z5 
< 




0.5 



2 4 6 8 10 12 14 
Radius (arcmin) 

EPIC MECS Ni Abundance profile 



1.5 



1 .0 



0.5 



□ PNS 
O PND 
★ MECS 



;> 



0.0 t 




2 4 6 

Radius (arcmin) 



Fig. 5. — Abundance profiles for Fe, on tlie top, and for Ni, on the bottom, obtained witfi 
tfie various EPIC cameras compared witfi tfie abundance profiles obtained with the MECS 
instrument on board BeppoSAX . Uncertainties are at the 68% level for one interesting 
parameter (Ax^ = 1). As in Fig. 3, crosses refers to MOSl, open triangles to M0S2, open 
squares to PN singles, open diamonds to PN doubles and full stars to MECS. 
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Fig. 6. — PN singles data for the 1' — 2' bin in the 6-8 keV band and the corresponding fit 
with an APEC model version 1.3.0, together with the corresponding ratios of data respect 
to models. 
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Fig. 7. — Temperature profiles obtained in the entire 0.5-10 keV band using one temperature 
(crosses) and two temperature models (open triangles, both the temperatures of the hot and 
cold component are shown) and, for comparison, the one obtained with the best fit model 
in the hard band (full circles), for PN singles data, on the top, and for M0S2 data, on the 
bottom. Uncertainties are at the 68% level for one interesting parameter (Ax^ = 1). 
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Fig. 8. — Ratio of the normalizations of the two temperature components, obtained by a 
two temperature model fitted in the 0.5-10 keV band, for PN single data (full circles) and 
M0S2 data (open triangles) . Uncertainties are at the 68% level for one interesting parameter 
(Ax^ = 1). 
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Fig. 9. — Fe abundance profiles obtained in the entire 0.5-10 keV band using one temperature 
(crosses) and two temperature models (open triangles) and, for comparison, the one obtained 
with the best fit model in the hard band (full circles), for PN singles data, on the top, and for 
M0S2 data, on the bottom. Uncertainties are at the 68% level for one interesting parameter 
(Ax' = 1). 
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Fig. 10. — Ni abundance profiles obtained in the entire 0.5-10 keV band using one tem- 
perature (crosses) and two temperature models (open triangles) and, for comparison, the 
one obtained with the best fit model in the hard band (full circles), for PN singles data, 
on the top, and for M0S2 data, on the bottom. Uncertainties are at the 68% level for one 
interesting parameter (Ax^ — 1). 



Table 1. Parameters values for one and two temperatures models in different energy 
bands for the inner rings of the Perseus cluster using the single events of the PN camera. 
All errors quoted are at the 68% level for one interesting parameter (Ax^ = 1 ) 



Bin 


Mod-Band 


Nh 


kTh 


EMh 


kTc 




Fc 


Ni 


xVd-o.f 






iiH-pow/D ^o-io.o Kev ) 




q Q1 +0.11 

"'•^■'--0.04 


A QQ+0-03 

4-99-0.09 






U.b2_o.oi 


1 rq+0.27 
i.OO_o.24 


QAK /Qnn 
o4o/ouy 






IT' in 

1 1 ^o-iU Kev J 




q QQ + 0.05 


A Q1+0.02 

4-91-0.03 






0.b2_o.oi 


1 .^0+0.27 
1.48_o24 


on A /TOtCi 

oo4/ <yy 






i 1 (U.o-iU kev j 


1 no+ooi 
1-09-0.01 


9Q+0.01 

O.Z»_g 01 


r 7C+0.01 

^- ' ^-0.01 






0-80-0.01 


l-^f5_0.15 


onQV / 1 onn 

zUo / / izyy 






Z L (^u.o-iu Kev ) 


1 1 c-l-O.Ol 

J--io_o.oi 


r nn+0-08 
5.00_o 01 


q r,7+0.02 
•^•"'-0.03 


9 n-i+o o^ 

Z.UO_o 01 


9 Qn+0.04 

2.yo_o.o7 


0-d9_o.oi 


1-Od_o.io 


1 7Q/t / 1 007 

1 / o4/ izy / 


1'. 


-2' 


1 1 H-pow/ D ^o-lo.o Kev j 




A ni+o.04 


1 1 q + 0.1 






n Rt^+O-Oi 
0.65_o.oi 


1 90+0.16 

J--^'->-0.16 


1114/ lUol 






IT / Q in 

i 1 ^o-iU Ke V j 




A 1 q + 0.06 


1 1 9+0.1 
11-^-0.1 






n R/i+0 01 
0.b4_o.oi 


1 9/1+0.15 
l-24_o.i7 


1 nni / 1 non 
luyi/ lUZU 






i 1 (^u.o-iu Kev ) 


1 19-1-0.01 


co+0.01 
0.00_Q 01 


19 1Q+O.OI 
1^-19-0.02 






0.«0_o.oi 


1 /19+O.IO 
l-4^_0.09 


OQ.QA /I KOn 

zoy4/ iozu 






zi ^u.o-iU Kev j 


1 17-1-0.01 
l-l'-o.oi 


A Qr;+0 02 
4-y5_o.03 


7 Q7+0.03 
'•"'-0.05 


^•U4-0.02 


A 77+0.08 

4-77_o.o3 


n 71 +0.01 


1 1R+O.O6 
1-lD-O.lO 


OQ71 /I c;i c 
Zo / 1/ lolo 


2'. 




1 rp 1 /I / o IOC 1 ATX 

il+pow/b (o-lo.o keVj 




A SR+0.07 


f^-^9-0.13 






n «i +0.01 
0-ol_o.oi 


n OC+0.15 

0-85_o.i7 


lL)/y/lL)Dci 






IT f'i in V/AA7'^ 
i i t^o-iU Ke V j 




4-99-0.04 


o.,+0.07 






A RI +0.01 

0-bl-o.oi 


n 01 +0.16 
0-91_o.i7 


1 n^n / 1 non 
lUoy/ lUZU 






IT fr\ K in i^/-.\7'\ 
i 1 (^U.o-iU Kev ) 


1 1^-1-0.01 


A 90 + 0.01 

4-^8-0.02 


Q 7^+0.01 






n 79+0.01 
0.72_o.oi 


1 9/ +0.13 

1.24_o^2 


OQ/IO /I Pii 

zo4y/ ioiy 






OT f'n rr; 1 n 

zi i^u.o-lU Kev j 


1 10-1-0.01 

i-i»_0.01 


r QC+0.04 

^•^^-0.03 


7 99+0.03 
' -^^-0.06 


9 11+0.02 

^•J-J--0.01 


9 oe+0.08 
^•»D_o.o2 


n R7+0.01 
0-b7_o.oi 


l-08_o.i2 


1 n7C / 1 1 7 

ly / 0/ loi / 


3' 


-4' 


IT+pow/b (3-13.5 keV) 




r (-9 + 0.08 

— u.u / 


6-44+°°4 

— U.U4 






0.56+°°i 

— U.Ui 


0-80+°1^9 


1061/1080 






IT (3-10 keV) 




'-'•^<5-0.09 


""J^-O.OS 






56+0 01 
^-"D-0.01 


0-87l°f, 


1062/1027 






IT-Fpow/b (0.5-10 keV) 


1 1 c+O.Ol 
-•^■^^-0.01 


4 qq+0.02 

^•^^-0.03 


6 93+0 01 
"•"^-0.02 






o-6itroi 


l-09l°:1^5 


2040/1527 






2T+pow/b (0.5-10 kcV) 


1 1Q-K0.01 


c C4+O.O2 


9C+0.03 


1 f;4+0.03 
i-O*-0.01 


n yq+o-oi 
u./y_o.o4 


0.5610.0? 


1 06+0.15 

J^-""-0.12 


1914/1525 


4'- 


-5' 


IT-Fpow/b (3-13.5 keV) 




c Q7+O.O8 
-0.09 


r f.7+0.08 

"3-0' -0.06 






1^2+0-01 

"-"J^-O.Ol 




1210/1090 






IT-Fpow/b (0.5-10 keV) 


1 10+0.02 

■'-•-'■"^-0.01 


r r7+0.02 


5 qq+0 01 
^•""-0.02 






"S5+0-01 

"-33-0.01 


o.8il°:l« 


1960/1521 






2T-hpow/b (0.5-10 keV) 


1 iq+o.oi 
i-iy_o.oi 




5.49tro? 


1 f.1+0.03 


60+0 01 

U.DU_o.o3 


^,s+o o3 
"-^•J-o.oi 


0.78l"„:l^ 


1877/1519 


5'- 


-6' 


IT-hpow/b (3-13.5 kcV) 






4.69tro^ 






0.4910-01 


^n+0.25 

U.OU_o 24 


1095/1072 






IT-Fpow/b (0.5-10 keV) 


1 10+0.01 
-'^•-'-"^-0.01 


c Gtc+0.02 

0.00_g 05 








0.5010.01 


1 14+0-21 
J.-J.^-0.18 


2004/1497 






2T-hpow/b (0.5-10 keV) 


1 10+0.01 

i-i»_0.01 


6.36tro^ 


/I RO+0.02 


1 07+0.03 


42+0 01 

U.*Z_0.03 


47+0.02 
u-*' -0.01 


1 24+0.14 

J--^*-0.23 


1898/1495 


6'. 


-8' 


IT-hpow/b (3-13.5 keV) 




6.78^°:1? 


7.44t"o;^J 






47+0.01 

"-^'-0.01 


12+0.23 

"-1^^-0.11 


1249/1289 






IT-hpow/b (0.5-10 keV) 


1 -ic+O.Ol 

^•-■^O-o.oi 


r, OC+0.03 
'J-°^-0.04 


8.04to-oi 






4q+o.03 
"-^"-0.01 


0.93l°?° 


2274/1637 






2T-|-pow/b (0.5-10 keV) 


1 91+0.01 

-^■^^-0.01 




7.5010-02 


1 90+0.01 
±.z,»_o 03 


63+0 02 

U-O'J-O.02 


47+0.01 

"-^' -0.01 


1 97+0.12 

.^-^' -0.23 


2009/1635 



Table 2. Parameters values for one and two temperatures models in different energy 
bands for the inner rings of the Perseus cluster using the the M0S2 camera. All errors 
quoted are at the 68% level for one interesting parameter (Ax^ = 1 ) 



Bin 


Mod- Band 




Ntt 

^ rl 


kTi, 


EMi, 


kT^ 




Fe 


Ni 


y2 /d f 


0.5'-!' 


IT+pow/b (3-12 keV) 






3.84iro^ 


5-36t°:°l 






0.64100? 


0.76+0-3^ 


385/323 




IT (3-10 keV) 






3.861°°^ 


r 07+0.02 

'J-'J ' _o 03 






0.64tn? 


0.8310-3° 


372/318 




IT (0.5-10 keV) 


1 


1 7-1-0.01 
1 —0.01 


00+0.01 

0.00_g Q2 


K qi+o oi 






O.Sltni 


l-04t°oJ^ 


810/479 




2T (0.5-10 keV) 


1 


9n-i-o.o2 

••^'J—O.Ol 


c -10+0.05 

d.lO_Q Qg 


9 09+0.03 

Z.OZ_g Qg 


9 00+0.02 


9C+0.25 
<J--^"_0.08 


0.74tni 


0.88^01^ 


fifi9/477 


l'-2' 


IT-hpow/b (3-12 keV) 






4 14+0.05 

—0 06 


12.2l°i 






0.72^001 


0.77+g-20 


43fi/393 




IT (3-10 keV) 






4.20iroi 








O.Tltni 


0.80tOQ;2i 


435/379 




IT (0.5-10 keV) 


1 


-lo-l-O.Ol 
•J-»_0 01 


3.67trol 


13.00^^2 






o.setHl 


l-09t°oJ9 


1228/540 




2T (0.5-10 keV) 


1 


-iq-FO.oi 
• j^y_o.oi 




4.09tH4 


2-9ltrol 


9.i5tn5 


o.83t"o;jl 


o.98i;3;?,J 


9fifi/538 


2'-3' 


IT-hpow/b (3-12 keV) 






4.96;ni 


io.ooi°:J^ 






0.64to°l 


0.49+J5-2? 


40fi/402 




IT (3-10 keV) 






r 1 7+0.05 
O.LI _g Qg 


9.83^007 






0.fi5t°Q°i 


0.59^0-^0 


415/383 




IT (0.5-10 keV) 


1 


1 7-1-0.01 
•J- ' _o 01 


4.35tro? 


10.96in? 






L*- ' ^ — 01 


i.J^'-' — O 11 


1106/544 




2T (0.5-10 keV) 


1 


91-|-0!01 

•^■■■-o.oi 


5 7S+0.05 

^- ' °-0.03 




9 09+0. 03 
^•''^-0.02 


40+0.08 


fiq+" "i 
^•"■^-0.01 


0.92tO;- 


841/542 


3'-4' 


IT-hpow/b (3-12 keV) 






r 70 + 0.09 

^•'°-0.10 


7 c;4+0.09 
' •O'i-0.07 






"•"J^-0.01 


o.8it°o;2^ 


407/406 




IT (3-10 keV) 






fi 1 5+°°^ 

O'-'-O-O.ll 


7 07+0.06 
' -"J ' -0.05 






^•"^-0.02 


1 03+0-2*5 

l.UO_Q 25 


415/380 




IT-l-pow/b (0.5-10 keV) 


1 


1 7+0.02 
-0.01 


r 9r+0.04 

^•^^-0.02 


7 qq+o ol 

' •^'-'-0.01 








-1 -17+0.22 
J^--"^ ' -0.15 


782/542 




2T+pow/b (0.5-10 kcV) 


1 


91+0.01 

■^■■■-o.oi 


"••J'J-o.io 


t- q-1 +0.03 
^•"JJ^-O.OS 


2.74tnf 


-1 47+0.14 

■■■•^'-0.06 




0.98tO;}« 


687/540 


4'-5' 


IT-hpow/b (3-12 keV) 






fi u+"-i2 


6.09t°;0« 






U.OO_Q Q]^ 


0.38t0;33 


412/401 




IT-l-pow/b (0.5-10 keV) 


1 


9R+0.01 
•^"-0.02 


5 fin+° °-'^ 

0-OU_o.04 


6.37t°;0^ 






t^8+o o^ 

^•3°-0.01 


1 29+0-21 

■■-•^■^-0.23 


701/538 




2T+pow/b (0.5-10 kcV) 


1 


97+0.02 
•^'-0.01 


fi.eelHi 


O-"J^_0.07 


9 qc+0.06 


-1 -10+0.12 

^^•^^'^-0.04 


n co+0.01 

U-JS-0.02 


-, 90+0.19 
^^•^•J-0.24 


666/536 


5'-6' 


IT+pow/b (3-12 kcV) 






-0.20 


^ no+0.06 

4-8J-0.05 






^8+o-03 

"■°°-0.03 


1 00+0-40 

l.UU_Q 42 


363/382 




IT-hpow/b (0.5-10 keV) 


1 


1 R+0.01 

•■LO-O.Ol 


fi 17+005 


5 01+° °2 

0-Ui_o.oi 






50+0 01 

u.oy_Q Q3 


1 fS8+0-28 
i.OO_Q 33 


650/518 




2T-|-pow/b (0.5-10 keV) 


1 


-ie+0.02 
•■'■''-0.02 


7 -1-1+0.07 
'•J^J^-0.19 


4 99+0.05 
*-'''_0.06 


0-1+0.15 
"^•"JJ^-O.IS 


o.89t°:J^g 


u.oy_o Q3 


1.0D_o 27 


639/516 


6'-8' 


IT-hpow/b (3-12 keV) 






' -^^-0.21 


7 on+0.09 

' -^^-0.09 






"■0^-0.03 


4fi+0-40 

'^•+"-0.39 


447/475 




IT-hpow/b (0.5-10 keV) 


1, 


90+0.02 
■^"-0.01 


fi ^n+0-06 


97+0.01 
' -0.02 






U.OO_Q QJ 


-I -17+0.28 
-■-•^ ' -0.26 


848/572 




2T-hpow/b (0.5-10 keV) 


1 


90+0.02 
•^"^-0.02 


7 qe+0.09 

' •yo_o.i4 




9 74+0.07 


i-46t°Q;^^3 


fifl+o oi 

U.OO_Q Q^ 


1 07+0-25 
-'-•'-"-0.26 


798/570 



